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The ground state phase diagram of the ID Hubbard chain with pair-hopping interaction is studied. 
The analysis of the model is performed using the continuum-limit field theory approach and exact 
diagonalization studies. At half-filling the phase diagram is shown to consist of two superconducting 
states with Cooper pair center-of-mass momentum Q = (BCS-770 phase) and Q — n (ry^-phase) and 
four insulating phases corresponding to the Mott antiferromagnet, the Peierls dimerized phase, the 
charge-density-wave (CDW) insulator as well as an unconventional insulating phase characterized by 
the coexistence of a CDW and a bond-located staggered magnetization. Away from half-filling the 
phase diagram consists of the superconducting BCS-770 and rj^ phases and the metallic Luttinger- 
liquid phase. The BCS-770 phase exhibits smooth crossover from a weak-coupling BCS type to a 
strong coupling local-pair regime. The rfr phase shows properties of the doublon (zero size Cooper 
pair) superconductor with Cooper pair center-of-mass momentum Q — n. The transition into the 
r^-paired state corresponds to an abrupt change in the groundstate structure. After the transition 
the conduction band is completely destroyed and a new r^-pair band corresponding to the strongly 
correlated doublon motion is created. 
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I. INTRODUCTION 

The problem of a crossover from a weak-coupling BCS 
picture of Cooper pair formation 0] to a Bose-Einstein 
(BE) condensation of preformed local pairs has drawn 
special interest during the last two decades (see for a re- 
view Refs. ). Increased renewed interest to this prob- 
lem mainly comes from experimental observations re- 
garding the unusual properties of high-T c materials. Par- 
ticularly important in this respect is the extreme short (of 
the order of one lattice spacing) coherence length of the 
pairs in the superconducting state [Q and a pseudo-gap 
structure in the normal-state density of states of under- 
doped cuprates ||. 

Leggett H was the first to point out that BCS super- 
conductivity can be continuously connected to BE con- 
densation by increasing the two-particle attraction be- 
tween electrons from weak to strong coupling. The size 
of Cooper pairs shrinks continuously until spatially well 
separated bosons form, which undergo Bose condensation 
at a sufficiently low temperature. Nozieres and Schmitt- 
Rink have extended Leggett 's work to lattice electrons 
and analyzed the crossover between the BCS and the BE 
limit as a function of coupling strength. The attractive 
(U < 0) Hubbard model 



H 



*5>Lc„ +lj(T + h.c.) + ^ Vp T (n)^(n) (1) 



was usually considered to describe the evolution from the 
BCS type pairing to the local pair (composite boson) 

rt 



limit lq-17 . Here a {c n a ) is the creation (annihila- 



tion) operator for an electron with spin a at site n and 
p a {n) = c+ c . 



A rather different realization of local-pairing is the 77 
-pairing mechanism of superconductivity introduced by 
Yang |L8| . Yang discovered a class of eigenstates of the 
Hubbard Hamiltonian which have t he p roperty of off- 
diagonal long-range order (ODLRO) |19[] , which in turn 
implies the Meissner effect and flux quantization p0| , pl| . 
These eigenstates are constructed in terms of doublon 
(on-site singlet pair) creation operators and can be shown 
to be of very high energy compared to the global ground 
state 18 1 . Following Yangs notations we call these 77- 
paired states. On a bipartite lattice one can consider two 
different realizations of the 77-paired state, constructed in 
terms of zero size Cooper pairs with center-of-mass mo- 
mentum equal to zero (770-pairing) and 7r (^-pairing), 
respectively. Yang also proved that these states can- 
not be ground states for the Hubbard model with fi- 
nite interaction |l8[ |. The 770-superconductivity is real- 
ized in the Hubbard model only at infinite on-site at- 
traction (22). Shortly after Yang's paper, Efiler, Kore- 
pin, and Schoutens proposed the supersymmetric exten- 
sion of the Hubbard model showing a true ODLRO and 
770-superconductivity in the ground state for a finite on- 
site interaction p3|. Later on several other integrable 
generalizations of the Hubbard model, with 770 type or- 
dering in the ground state were proposed The im- 
portant common feature for all integrable models of 77- 
superconductivity is the high SU(2) eg) SU(2) symmetry 
which is ensured by the strongly correlated kinematics of 
electrons on a lattice. 

Another model with a kinematical mechanism for the 
formation of Cooper pairs is the Penson-Kolb (PK) 
model |25|. The Hamiltonian of the PK model contains, 
in addition to the usual one-electron hopping term, a 



1 



term that moves singlet pairs of electrons from site to 
site and in the one-dimensional case is given by 



ft = "^(iA+l.a + h - c -) 
n, a 

+ W^{q ] {n)q i {n) + h.c). 



(2) 



where q a (n) 



= A 



In the case W < the PK 



model (^) describes a continuous evolution of the usual 
BCS type superconducting state at \W\ <C t into a local 
pair ?7o-type state at \W\ ^> t p6| ]. More important is 
that in the case of "repulsive" pair-hopping interaction 
(W > 0), the ^-paired state is realized in the ground 
state of the PK model for W > W c ~ 2t 

It is notable that the pair-hopping term and the Hub- 
bard term could be obtained from the same general tight- 
binding Hamiltonian (30| by focusing on selected terms of 
the two-particle interaction. Indeed the same matrix ele- 
ment of the electron-electron interaction potential V(r), 
V(n, m, k, I) — (n, m\V(r)\k, I) which gives rise to the on- 
site Hubbard interaction for n — m = k = I leads to the 
pair-hopping amplitude W for n — m and k = I = n± 1. 
Although originating from the same two-body potential 
the Hubbard and the pair-hopping couplings represent 
different types of correlations in the electron system. If 
the Hubbard term describes on-site correlations, the site- 
off-diagonal pair-hopping term describes part of the so- 
called "bond-charge" interaction. The sign of Coulomb- 
driven on-site and pair-hopping interaction is typically 
repulsive U, W > 0. However, we will treat the parame- 
ters U, W as the effective ones, assuming that they in- 
clude all the possible renormalizations. In particular, 
contribution from the strong electron-phonon coupling 
or from the coupling between electrons and other elec- 
tronic subsystem could give the effective attractive on- 
site (U < 0) interaction j§]. The "attractive" (W < 0) 
pair-hopping term can originate from the coupling of elec- 
trons with intermolecular vibrations |3l| ] , or from the on- 
site hybridization term in a generalized periodic Ander- 
son model p^ |. 

In this paper we investigate the interplay between 
these different sources for local pair formation: the site- 
diagonal Hubbard attraction and the site-off-diagonal 
pair-hopping interaction. We consider the extended Hub- 
bard chain with pair-hopping interaction i.e. the so 
called Penson-Kolb-Hubbard (PKH) model (3^]. In one- 
dimension the Hamiltonian reads: 

n,a n 

+ Wj2(^( n )Ql( n ) + h.c.)+^P<r{n). (3) 

n n.a 

There are N e particles, Nq sites and the band filling v — 
N e /2No is controlled by the chemical potential [i. In the 
absence of the W term the Hamiltonian (||) corresponds 



to the Hubbard model, while in the absence of the U 
term it reduces to the Penson-Kolb model. 

The PKH model has been investigated mainly in the 
36|| . Recently Robaszkiewicz and 
< PKH model by means of the 



34 



case U > |3§| 
Bulka studied the U 
Hartree-Fock approximation and the slave-boson mean- 
field method j37|. In this paper we focus on the na- 
ture of phase transitions in the ground state of the PKH 
model between the following local-pair ordered phases in 
p7| , ^6| : the ?7o-superconductor, the CDW phase, and the 
TyTr-superconductor. 

The paper is organized as follows: In Sect. II we review 
the model and its symmetries. We study the ground state 
phase diagram of the PKH chain using the continuum- 
limit bosonization approach (Sect III) and exact Lanczos 
diagonalization for chains up to L = 12 sites (Sect. IV). 
Finally, Sect. V is devoted to a discussion and to con- 
cluding remarks. 



II. REVIEW OF THE PKH MODEL 

There are two important aspects distinguishing the 
PKH model from the Hubbard model: the symmetry 
of these models and the nonlocal character of the pair- 
hopping interaction. 

Let us first consider the symmetry aspect. The PKH 
model and the Hubbard model are characterized by the 
same SU (2)-spm symmetry. The difference lies in the 
symmetry of the corresponding charge sectors - SU (2) 
in the case of the 1/2-filled Hubbard model and U(l) in 
the case of the PK model. This can be easily seen for a 
strong Hubbard attraction \U\ ^> \W\,t where a large gap 
of order \U\ exists in the single-particle excitation spec- 
trum. Projecting the system on the subspace excluding 
single occupancy of sites and using the standard second- 
order perturbation theory with respect to t 2 /\U\ |38| one 
obtains the following effective spin-I/2 XX Z spin chain 
Hamiltonian 

H = E a + JY$(Vo<n)Vo(n + 1) + h - c -) 



+ Ar, z (n)r] z (n + 1)}, 



where E = -±N \U\, 



4t 2 /\U\ 



At 2 

J = 2W ~WY A " |2VF-4t 2 /|[/|| 
and the pseudospin operators are 

r)+(n) = 4 )T 4 4 , »7o (n) = c„ a c„, T , 

^oW = (4,T c ».T+4a c n4- 1 )/ 2 - 



(4) 



(5) 



(6) 



As we see the charge sectors of the half-filled PKH 
and Hubbard models are governed by the U(l) (A ^ I ) 
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and the SU(2) (A = 1) symmetry of the equivalent 
Heisenberg model respectively. In the case of the half- 
filled Hubbard model, due to the Si7(2)-charge symme- 
try, the superconducting order mixes with the charge 
density order. The singlet superconducting (SS) and the 
charge density wave (CDW) correlations show an identi- 
cal power-law decay in the infrared (large distance) limit 

s 

For W (A ^ 1) the spin rj = 1/2 antiferromagnet, 
with Hamiltonian (Q) is in a gapless planar XY phase for 
A < 1 and in the gapped Neel phase for A > 1 . The 
Neel phase is realized for < W < At 2 /\U\ and corre- 
sponds to the insulating long-range-ordered (LRO) CDW 
state in terms of the initial electron system. The gapless 
planar XY phase is realized for W < and W > At 2 /\U\. 
The tendency towards in-plane magnetic ordering reflects 
to the superconducting ordering within the initial elec- 
tron system. 

In the ultimate limit t 2 /\U\ -> the PKH model is 
equivalent to the spin- 1/2 XY model. Electrons only 
appear in singlet pairs on the same site, the interaction 
simply acts as a hopping term for these pairs. The model 
is ry-superconducting by construction. The CDW corre- 
lations 



Al 



(ii) 



(7) 



decay faster then the ^-superconducting correlations 

(4{Q)^{n))~n- l '\ (8) 

This picture holds true for arbitrary W . However, if 
the jyo-ordering is realized for W < 0, then in the case 
of a "repulsive" (W > 0) pair-hopping coupling , the 
^-superconducting correlations 



(9) 



have to dominate. In deriving (J9J) the unitary transforma- 
tion 77(j~(n) — > (—l) n rjQ (n) = r]+(n) has been used which 
changes the sign of the transverse exchange Jj_ — » — J± . 

Although the XY model is invariant with respect to 
a sign change of the coupling constant, W — ► — W is 
not a symmetry of the PKH model. Therefore the way, 
the system approaches its limiting behaviour at \W\ ^> 
t 2 /\U\ is genuinely different for negative and positive W. 

The origin of this difference arises from the site- off- 
diagonal nature of the pair-hopping term. As far as the 
essence of the W <-> —W asymetry is connected to two 
different possibilities for Bose condensation of Cooper 
pairs with momentum Q — and Q = n it is conve- 
nient to rewrite the Hamiltonian (^) in momentum space 
as 

H = -2t£ C L%* cos ( fc ) + E V (Q) A Q A Q- ( 10 ) 

fe.CT Q 

Here V(Q) = U + 2Wcos(Q) and 



is the creation operator for a pair of electrons with op- 
posite spins and total momentum Q. 

In the weak-coupling (large bandwidth) limit U, W <C t 
only scattering processes involving states near the two 
Fermi points kp — ±ttv are relevant. Therefore, the scat- 
tering processes with Q ~ (forward scattering without 
and with spin flip) are characterized by the effective cou- 
pling constant Vfs(Q kp) ~ U + 2H^ and the umk- 
lapp scattering with Q ~ n by the effective amplitude 
coupling constant VjjmkiQ ~ t) — U — 2W. The umk- 
lapp scattering is irrelevant at noncommensurate band 
filling v 7^ 1/2. Therefore, at v ^ 1/2 the weak-coupling 
phase diagram of the PKH model has to be similar to 
that of the Hubbard model with an effective interaction 
U eff = U + 2W @-||. At half-filling the umklapp 
scattering is relevant in the case of repulsive interaction 
(U, W > 0) and leads to the CDW type ordering in the 
PKH model at < U < 2W (36). 

In the forthcoming section we will use the continuum- 
limit bosonization approach to study the ground state 
phase diagram of the PKH model. However, already 
the above presented qualitative analysis indicates the 
absence of the ^-superconducting phase in the weak- 
coupling phase diagram. 

The particularity of the model in the case of strong 
pair-hopping interaction (\W\ ^> t) could be observed 
even in the simple case of two particles with opposite 
spins on the lattice. As far as the total momentum of the 
system Q tot is conserved, one can treat each Q-sector of 
the Hilbert space independently. Since the ground state 
belongs to the sectors with Qtot = or Qtot = 7r, we 
restrict ourselves to these sectors only. The states with 
Qtot i= or 7r exhibit a broken time reversal symmetry 
and are excited states. 

Let us first consider the U = case. At W = the 
ground state energy is Eq = —At and the total momen- 
tum is Qtot — 0. The probability to find the on-site pair 
in the ground state is 1/L. Eigenstates in the Qtot — ^ 
sector correspond to the highly excited states with the 
energy E = 0. At |W|/t — > oo the r]o pair given by the 
wave function Aq\Q > has the energy E Vo = 2W while the 
r\n pair given by the wave function A\ |0 > has the energy 
E ri7r = —2W. Therefore, for W < 0, the groundstate al- 
ways remains in the Qtot = subspace and its energy con- 
tinuously varies from -At (W = 0) to -2\W\ (\W\ > t). 
The probability to find the on-site pair increases contin- 
uously up to 1 (|T^|/i — + oo). In contrast, for W > 0, the 
energy of an 770-paired state goes to E Vo = 2W > while 
the energy of an r^-paired state E Vjr = —2W < 0. Thus, 
with increasing W/t, at some critical value of the pair- 
hopping coupling W c > 0, the total momentum of the 
system in the ground state should change from Qtot — 
to Qtot — 0. In the Qtot = t sector the contribution 
of the one-particle hopping term to the energy vanishes. 
Below we will refer to this phenomenon as the collapse 
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of a one-particle band. On the other hand the contribu- 
tion of the derealization energy for an r)„ pair is equal 
to —2W. Therefore one can roughly estimate the critical 
value to be W c — —It. 

We now consider the effect of the Hubbard interaction. 
It is clear that the on-site repulsion U > increases the 
difference between the unpaired state in the Qtot = 
sector and the ^-paired state in the Qtot — ^ sector. 
Therefore the critical value has to increase in the case of 
on-site repulsion. Contrary, the attractive Hubbard in- 
teraction U < supports the on-site pairing, reduces the 
effect of the t-term and therefore the difference between 
the lowest energy states in the Qtot = and Qtot = t sec- 
tors. Therefore, in agreement with the strong-coupling 
analysis the critical value has to strongly decrease in the 
case of on-site attraction. As we show below, the picture 
remains qualitatively similar in the case of many particle 
systems. The essence of the transition into an ^-paired 
state consists of a collapse of the one-particle band and 
the creation of a strongly correlated two-particle 77-pair 
band. 



III. WEAK-COUPLING PHASE DIAGRAM. 

The field theory treatment of ID systems of corre- 
lated electrons is based on the weak-coupling approach 
(\U\, \ W\ <i). Assuming that the low energy physics is 
controlled by states near the Fermi points ±fcp (fefao = 
ttv, where ao is the lattice spacing) we linearize the 
spectrum around these points and obtain two species 
(for each spin projection a) of fcrmions, R a (n) and 
L (T (n), which describe excitations with dispersion rela- 
tions E(p) = ±vfP- Here, vf = 2taosin(iri/) is the Fermi 
velocity and the momentum p is measured from the two 
Fermi points. More explicitly, one decomposes the mo- 
mentum expansion for the initial lattice operators into 
two parts centered around ±fcp to obtain the mapping: 



l R a {n)+e-™ vn L a {n), 



(12) 



where the fields R a (n) and L a (n) describe right and left- 
moving particles, respectively, and are assumed to be 
smooth on the scale of the lattice spacing. This allows us 
to introduce the continuum fields R a (x) and L a (x) by 



Rain) 



a R a (x = na ), 
oqL^x = na ). 



(13) 



In terms of the continuum fields the free Hamiltonian 
reads: 



Ho 



-iVF^y] J dx[: R\d x R a : - : L\d x L a :] 



(14) 



which is recognized as the Hamiltonian of a free massless 
Dirac field. Here and below :...: denotes normal order- 
ing with respect to the ground state of the free electron 
system. 



The advantage of the linearization of the spectrum 
is twofold: the initial lattice problem is reformulated 
in terms of smooth continuum fields and - using the 
bosonization procedure - is mapped to the theory of two 
independent (in the weak-coupling limit) quantum Bosc 
fields describing charge and spin degrees of freedom, re- 
spectively. 

In terms of the continuum fields the initial lattice op- 
erators have the form: 

Pa{n) -> ao{(Ji?,a(aO + Jl,<t(x)) 

+ (e- i2 ™ n Rl{x)L a {x)+h.c.)}, (15) 

+ (e~ l ^ 2n+ ^Rl{x)L a {x) + h.c)y (16) 

where 

J R< a =: R ] a (x)R a {x) :, J L . a =: Ll(x)L r7 (x) : . 

The second step is to use the standard bosonization ex- 
pressions for fermionic bilinears jy] : 

-i^[:Rld x R„ : - : L\d x L a : ] 

a 

\{(d x 9 c f + (d x cf> c ) 2 } + \{{d x e s ) 2 + (d x <p s ) 2 }, (17) 



Jr,o + Jh,o 
Jll,a — Jh,a 

R\(x)L a {x) 



1 



V27T 
1 

—i 
2nao 



[{d x (/>c) + cr(d x (f> a )] , 

[(d x e c ) + o-(d x e s )] 

exp (zv / 27r((/) c + o-<f> s )). 



(18) 
(19) 
(20) 



We thereby obtain 

Pa{n) -> a \— == \{d x <j) c ) + a(d x (p s )] 



. 2?r 

■ sin (v / 27r( 



— 2nnv 



na 



>}■ 

q a {n) -> a \ —= cos(irv){(d x <j) c ) + cr{d x 4> s )} 
V 2tt 

:sm(irv)[(d x 6 c )+a(d x 6 s )] 

a(f> s ) - (2n + l)m>)\, 



(21) 



+ 



'2-k 
sin (\2~n( 



(22) 



Here scalar fields 4>c,s(x) describe the charge and the spin 
degrees of freedom and fields 9 c ,s(x) are their dual coun- 
terparts: d x 9 c ^ s = Tl c ,s where n cs is the momentum con- 
jugate to the field 4> c , s - Inserting the relations (|2l])-(|2^) 
into (0), by virtue of the smooth variation of the fields 
(j) c . s (x) on the scale of lattice spacing one can convert 
sums into integrals ao — > J dx and take the contin- 
uum limit. 



4 



All terms in the obtained continuum-limit Hamiltonian 
that contain the rapidly oscillating phase factors e ±l2n7rl/ 
and e ±l4n71 ' l/ will drop out except at v = 1/2, where terms 
containing e ±l4n7TL ' = l ; are not oscillating and have to be 
kept. Performing this procedure, after rescaling of fields 
and lengths, the continuum-limit version of the Hamilto- 
nian (^) acquires the following form: 



(23) 



where 



n s 



d.r 



2n 2 a 2 



{y[(^ s ) 2 + (^ s ) 2 ] 
(y/8nK a <p 3 (x))} 



(24) 



describes the spin degrees of freedom. 

The charge degrees of freedom for v ^ 1/2 are de- 
scribed by the free scalar field 



H c = ^Jdx [(d xlfc ) 2 + (d x ti c ) 2 ] , (25) 
and for v = 1/2 by the quantum sine-Gordon field 
H c = J dx{^ [(d x <p c ) 2 + d x $ c ) 2 ] 

+ cos WtorK c <p c (x))\. (26) 

All CLq ) 



Here we have defined 



Kr 



777. r 



K. 



U-2W 



1 



2irt 
1 



2" 

Z7 + 2VT 
2nt ' 



1 [/ 

2 _ 

1 C7 
2 _ 



2PT 



27rf 



2PT 



2?rf 



(27) 

(28) 

(29) 

(30) 
(31) 



The mapping of the Hamiltonian (||) into the quantum 
theory of two independent charge and spin Bose fields, 
allows to study the ground state phase diagram of the 
initial electron system, using the far-infrared properties 
of the bosonic Hamiltonians j2^)-(26). 



A. Renormalization group analysis 

Let us first consider the half-filled band case when both 
the spin and the charge sectors of the system are gov- 
erned by the quantum sine-Gordon (SG) fields and 
(p6|). The infrared behavior of the SG Hamiltonian TL C . S 
is described by the corresponding pair of renormalization 
group (RG) equations for the effective coupling constants 



M c ( s ) and K, 




2(K-1) 



FIG. 1. Renormalization-group flow diagram; the arrows 
denote the direction of flow with increasing length scale. The 
circles correspond to the starting points in the case of a gap- 
less (A) and a gapped (Ai) spin channel. The crosses corre- 
spond to the starting points for two different gapped sectors 
in the charge channel with {tp c ) = (t3i) and (tp c ) = ir/8K c 
(B 2 ). 



dM c(s) {L) 

dL 

dK c{s) (L) 

dL 



= -2{K c{s) {L)-l)M c{s) {L) 



= -~Ml s) (L), 



(32) 



and M c(s) (0) = 
2T) describes the 



where L = ln(a ), K c(s) (0) = 1 + ±g c( 
m c r s ). Each pair of RG equations ( | 
Kosterlitz-Thouless transition |Q in the charge and spin 
channels. The flow lines lie on the hyperbola 



4(#eM " I) 2 - M 2 C[S) = £ {s) = g 2 c(s) 



c(a)> 



(33) 



and - depending on the relation between the bare cou- 
pling constants g c [ s ) = 2{K c i s \ — 1) and 7n c ( s ) - exhibit 
two different regimes (see FigJl]): 

For g c i s \ > \m c f s \ | we are in the weak coupling regime; 



the effective mass M, 



0. The low energy (large 



distance) behavior of the gapless charge (spin) excitations 
is described by a free scalar field. The corresponding 
correlations show a power law decay 



/ i V2TrK*tp ( x ) — i V^TrlC'ip ( x ' ) \ 



\X — X 



-K* 



( e i^/2w/K'e(x) e -iy/2n/K*0(x')\ _ u _ x >\- l / K * 



(34) 
(35) 



c( S ) 



and the only parameter controlling the infrared behav- 
ior in the gapless regime is the fixed-point value of the 
effective coupling constants K*( a y 

For g c r s \ < |t77 c ( s ) | the system scales into a strong cou- 
pling regime: depending on the sign of the bare mass 
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m c (s), the effective mass M c ( s ) 



±00, which signals the 
crossover into a strong coupling regime and indicates the 
dynamical generation of a commensurability gap in the 



excitation spectrum. The field 
the vacuum expectation values | 



£c( s ) 



gets ordered with 



■ , o j? i m c(s) > 0) 

(ip) = { V 8/ M s ) 

(m c(s) < 0) 



(36) 



The ordering of these fields determines the symmetry 
properties of the possible ordered groundstates of the 
fermionic system. 

Using Eqs. ©-(fl]) and ([|| one easily finds that 
there is a gap in the spin excitation spectrum (M s — > 
— oo) for U + 2W < 0. In the spin-gap sector the ip s field 
gets ordered with vacuum expectation value ((p s ) = 0. In 
the sector U + 2W > the spin excitations are gapless 
(M s — > 0) and the low-energy properties of the spin sector 
are described by the free Bose field system with the fixed- 
point value of the parameter K* = 1. 

At half-filling the charge sector is gapless for U, W < 
and along the line U = 2W for U, W > 0. At U, W < 
the low-energy properties of the gapless charge sector are 
described by the free Bose field Hamiltonian (25) with the 
fixed-point value of the parameter 



K* 



1 + V2UW/2nt. 



(37) 



For U, W > the line U — 2W (m c = 0) corresponding to 
the unstable fixed-point line m c — 0, K c — 1 < (see Fig. 
|l|) . Here the infrared properties of the gapless charge sec- 
tor are described by the free massless Bose field with the 
bare value of the Luttinger liquid parameter K c . More- 
over, the line m c = 0, K c - 1 < (U, W > 0, U = 2W) 
separates two different insulating (charge gapped) sec- 
tors of the phase diagram: for U > max{2W^, 0}, m c < 0, 
M c — > — oo and therefore the ip c field gets ordered with 
vacuum expectation value (<p c ) = 0; for 2W > max{[/, 0} 
m c > M c — » +oo and therefore the ip c field gets ordered 
with vacuum expectation value (tp c ) = y/n/8K c . 

Away from half-filling the charge sector is gapless and 
is described by the free massless Bose field (p5|). The 
corresponding correlations (34 )-(^5|) show a power law 
decay at large distances with critical indices determined 
by the bare value of the coupling constant K c ( p7j ) . The 
spin channel remains massive at U + 2W < and gapless 
for U + 2W > 0. 



B. Correlation functions 

To clarify the symmetry properties of the ground states 
of the system in different sectors we introduce the follow- 
ing set of order parameters describing the short wave- 
length fluctuations of the site-located charge and spin 
density 



&CDw(n) = e l2 ™ n J2 a p a (n) ^ 

sin (y/2-KK c (p c (x)) cos ( y \/2 / nK s ip s {x)) , at v = l/2 
e W27rK cVc (x) cog (^/2lrfr st p s {x)) , at 



30S (^/2irK c (p c {x)) sin ( y ^/2iiK s (p s (x)') , at v = 1/2 
3 W2*k cVc (x) gin (^2irK~ s Lp s (x)) , at i^l/2 

and two superconducting order parameters correspond- 
ing to singlet (Ass) and triplet (Ays) superconductivity: 

A SS {x) = R\{x)L\{x) - r\(x)l\(x) - 

cxp ^i^j^-O c (x) \ cos (y/2irK s ip s (x)), 

A TS (x) = R\(x)l\(x) + r\(x)l\(x) 



(40) 




<'xp ( i\/ —6 c (x) ) sin (y/27rK s <p s (x)). 



(41) 



In the particular case of a half-filled band, we consider an 
additional pair of order parameters, corresponding to the 
short wavelength fluctuations of the bond-located charge 
and spin density 



A Dmer (n) = (-1)" ^(^(n) + # („)) 



COS (y / 2TTK c (f c (x)) COS (y / 2TTK s (f s (x)), 



A Bd -SDw(n) = (-l) n J2 a (&( n )+€(n)) 



sm 



(y/2ivK c (p c (x)) sin (\/2TrK s (p s (x)), 



(42) 



(43) 



With these results for the excitation spectrum and the 
behavior of the corresponding fields Eqs. (34)-([36"|) we 
now discuss the weak-coupling ground state phase dia- 
gram of the model (||) . 



C. The half-filled band case. 

At half-filling the weak-coupling phase diagram of the 
PKH model consists of the following sectors (see Fig. 2). 

A. The Singlet Superconducting sector: 
U < and W < 0; 
M S ^Q, (yj,) = 0; M c = 0, K*>1. 

This sector of the coupling constants corresponds to 
the strong-coupling regime in the spin-channel and gap- 
less charge channel. The dynamical generation of a spin 
gap accompanied by the ordering of the field ip s with 
vacuum expectation value (<^ s ) = leads to a complete 
suppression of the SDW, Bd-SDW, and TS instabilities. 
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FIG. 2. The groundstate phase diagram of the ID Pen- 
son-Kolb-Hubbard model at half-filling. The dot-dashed line 
marks a transition into the ^-superconducting phase. Solid 
lines separate different phases: l.Dimer - long range ordered 
(LRO) dimerized (Peierls) phase. 2. CDW - LRO charge 
density wave phase; 3. SDW + Dimer - insulating phase 
with an identical power-law decay of spin-density-wave and 
Peierls correlations. 4. CDW + Bd-SDW - insulating phase 
with an identical power-law decay of CDW and bond-located 
SDW correlations. 5. BCS - ^-superconductor - singlet su- 
perconducting phase exhibiting a continuous evolution with 
increasing \U\/t and/or IW^/f from the standard BCS limit 
into the local pair no state. 

The CDW, Dimer, and SS instabilities survive and show 
a power-law decay at large distances 

{A C Dw{x)A C Dw{x')) ~ (A Dimer 

(x)A 

Dimer 

CO) 

~\x-x'\- K "% (44) 

(Ass(x)Ass(x')) ~\x-x'\- 1/K ■ (45) 

where K* is the fixed-point value of the parameter K c . 
For U, W < 0, K* > 1 and the SS instability dominates 
in the groundstate. 

Note that, at W = or U = 0, K* = 1 and the CDW, 
Dimer, and SS correlations show an identical power law 
decay at large distances. At W — coexistence of the 
CDW and superconducting instabilities in the ground 
state reflects the high SU(2) <g) SU(2) symmetry of the 
1/2-filled Hubbard model p9| . In the case of the Penson- 
Kolb model (U — 0), due to the charge-JJ(l) symmetry 
of the pair-hopping term, there is no symmetry reason 
which guarantees scaling to the SU(2) invariant fixed 
point K* = 1. As it was shown by Affleck and Marston, 
the higher order RG corrections give K* > 1 indicat- 
ing the dominating character of the SS instability in the 



ground state of the 1/2-filled PK model for arbitrary 
W < @. 

B. The CDW sector: 

U < 0, < W < -U/2. 

M s £ 0, (p a ) = 0; M c ? 0, (<p c ) = ^JW C . 

For W > the commensurability gap in the charge 
degrees of freedom opens. For < W < —U/2 the spin 
sector remains massive. Ordering of the field tp s with 
vacuum expectation value (tp a ) = leads to a complete 
suppression of the SDW, Bd-SDW, and TS instabilities. 
Ordering of the field ip c with vacuum expectation value 
((f c ) = y/Tv/SKc leads to a suppression of the SS and 
Dimer correlations. The CDW correlations show a true 
long-range order 

(A C dw{x)A C dw{x')) — const (46) 

in the ground state. Therefore in this sector of the phase 
diagram the system displays the properties of the CDW 
insulator. 

C. The Peierls (dimerized) sector: 
U > 0, < U < -2W. 

M s £ 0, {<p a ) = 0; M c £ 0, ( Vc ) = 0. 

Ordering of the field ip s with vacuum expectation value 
(ip s ) = leads to a complete suppression of the SDW, Bd- 
SDW, and TS instabilities. Ordering of the field <p c with 
vacuum expectation value (ip c ) = leads to a suppression 
of the SS and CDW correlations. The Dimer correlations 
show a true long-range order 

(A Dimer (x)A Dimer (x')) - const (47) 

in the ground state. Therefore in this sector of the phase 
diagram the system is a dimerized (Peierls) insulator. 

D. The (CDW + Bd-SDW) sector: 
W > 0, -2VF < U < 2W. 

m s = o,k; = i, m c ^ 0, (^ c > - ^JW C . 

The generation of a gap in the charge excitation spec- 
trum, accompanied by the ordering of the field (p c with 
vacuum expectation value (ip c ) — ^tt/8K c leads to a 
suppression of the superconducting, SDW, and Dimer 
ordering. The CDW and Bd-SDW correlations show a 
power-law decay at large distances 

(Acdw(x)Acdw(x')) ~ (A B d-SDw{x)A B d-SDw{x')) 
-Ix-x'f 1 . (48) 

Therefore this sector of the phase diagram corresponds to 
the insulating phase with coexisting CDW and Bd-SDW 
instabilities. 

E. The (SDW + Dimer) sector: 
U > 2\W\. 

Ms = o,k; = i, m c 0, (cp c ) = 0. 
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FIG. 3. The ground state phase diagram of the ID Pen- 
son- Kolb-Hubbard model for v ^ 1/2. The dot-dashed line 
marks the transition into the ^-superconducting phase. The 
solid line separates the superconducting phase from the Lut- 
tinger metal phase characterized by an identical power-law 
decay of the SDW, CDW, Bd-SDW, and Dimer correlations. 

The generation of a gap in the charge excitation spec- 
trum, accompanied by the ordering of the field <p c with 
vacuum expectation value (ip c ) = leads to a suppres- 
sion of the superconducting, CDW, and Bd-SDW corre- 
lations. The SDW and Dimer correlations show a power- 
law decay at large distances 

~\x-x'\~ l . (49) 

Therefore this sector of the phase diagram corresponds 
to the antiferromagnetic insulating phase with coexisting 
SDW and Pcicrls instabilities. 



1. The non- half-filled band case. 

Away from half-filling and for U < the weak-coupling 
phase diagram of the PKH model consists of the following 
two sectors. 

Al. The SS sector: 
U + 2W < 0. 

M s f 0, (<p a ) = 0; M c = 0,K C > 1. 

In this sector of coupling constants the SDW, Bd- 
SDW, and TS instabilities are suppressed. The CDW, 
Dimer, and SS instabilities show a power-law decay at 
large distances 

(x) A Dimer 

~\ x -x'\~ Kc , (50) 



(A ss (x)Ass(x')) ~ \x-x'\- 1/Ke . (51) 

As it follows from @, for U = 2W < 0, K c > 1 and the 
SS instability dominates in the groundstate. 

Bl. The metallic Luttinger-liquid sector: 
U + 2W> 0. 

M s = 0, K* = 1, M c = 0,K C < 1. 

Both the charge and the spin excitations are gapless. 
All correlations show a power-law decay in the infrared 
limit. However, as far as K c < 1, the superconducting 
correlations 

(A ss (x)Ass(x')) (A TS (x)A TS (x')) 



decay faster then the density-density correlations 

(Ai(x)Ai(x')) ~\x — x'Y X ~ Ka , (53) 

where i = CDW,SDW,Dimer,Bd-SDW. Therefore in this 
sector the ground state of the PKH model shows proper- 
ties of the Luttinger liquid phase with a weakly dominat- 
ing tendency towards the density-density type ordering. 

To summarize this section, we have presented the 
weak-coupling groundstate phase diagram for onc- 
dimensional PKH model for arbitrary U and W 
(\U\, | W\ -C t). We have shown that the model has a very 
rich phase diagram including at half-filling, the singlet- 
superconducting phase (U, W < 0) and four different in- 
sulating phases corresponding to the Mott antiferromag- 
net (U > 2\W\), the Peierls dimerized insulator (0 < U < 
-2W) and the CDW insulator (W > 0, U < -2W) and 
an unconventional insulating phase characterized by the 
coexistence of the CDW and the bond-located staggered 
magnetization Bd-SDW (W > 0, U < 2\W\). The possi- 
bility of bond-located ordering results from the site-off- 
diagonal nature of the pair-hopping term and is a special 
feature of the half-filled band case. Away from half-filling 
the phase diagram consists of the singlet-superconducting 
phase (U + 2W < 0) and the metallic Luttinger-liquid 
phase (U + 2W > 0). 

Note the absence of the ^-superconducting phase in 
the weak-coupling phase diagram. This results from the 
finite-band nature of the transition into an f^-paired 
phase J2?j] . Therefore in the used infinite-band bosoniza- 
tion approach, such a transition could not be traced. In 
the forthcoming section we will use exact Lanczos diago- 
nalizations for chains up to L = 12 sites to study in detail 
the nature of the transition into the Ty^-superconducting 
phase. 

IV. NUMERICAL STUDIES 

In this section we study the nature of the transition 
into the ^^-superconducting phase using exact diagonal- 
izations for chains up to L = 12 sites. 
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A. The half-filled band case 



We start with the W/t dependence of the groundstate 
energy for different U < 0. 

In Fig. ||a we have plotted the groundstate energy 
(Eos) shifted by ±N e \U\ of the L — 10 chain at half- 
filling as a function of W/t, for selected values of U/t. 
We clearly observe a different behavior of the groundstate 
energy for W > and W < 0. For W < the ground 
state energy continuously approaches the XY model in 
the limit | W \— * oo. There is no trace of any additional 
transition for W < 0, in agreement with the density- 
matrix renormalization group (DMRG) results for U = 
p9[ . However, for W > 0, the groundstate energy shows 
a non-monotonic behaviour. We observe that above some 
critical value of the pair-hopping coupling W c the en- 
ergy is already almost linear in W and approaches the 
groundstate energy of the XY model. Hence for W > W c 
the one particle hopping term is almost frozen out. This 
change of behaviour is attributed to the transition to the 
ry^-superconducting state p7[ . The accurate definition 
for W c will be given below; however, the numerical data 
presented in Fig. [|a already clearly indicate the renor- 
malization of the critical value of the pair-hopping cou- 
pling W c by the on-site Hubbard attraction. At U — 
W c ~ 1.8* @ and it reduces to W c ~ OAt at U/t = -8. 

Another way to visualize this transition is to show the 
pairing phenomenon in real space. The corresponding 
quantity is the expectation value to find an on-site pair 
and is given by 



1 L 

Prs = ~£ / A/VrA-q)- 



(54) 



P rs is shown in Fig. ^b. For W < the pairs appear 
continuously in the system. In the opposite case W > 0, 
W <C t and \U\ < At the tendency to pairing is even 
reduced. This "effective repulsive" character is absent 
for U < —it, but remains up to W < 0.8t for U = 
and up to W < OAt for U = —It. Wc attribute this 
behaviour to a finite size effect where the energy scale 
coming from the exponentially small charge gap in W/t 
is not traced in the L = 10 chain. For U < —it the 
tendency to pairing is present for W > 0. As is seen 
from Fig. |]b for arbitrary U < 0, the almost fully paired 
state is realized in the ground state of the half-filled PKH 
model for a finite value of the pair-hopping amplitude 
W > W c . However, it should be stressed that due to 
the quasi-bosonic character of the pairs, the weight of 
unpaired particles remains extremely small but finite even 
for W ^ W c - The critical value of the pair- hopping 
amplitude corresponding to the transition into an rj^- 
superconducting state W c is straightforward to find when 
the number of electrons is N e = 2(2n + 1). The level 
crossing phenomenon is generic in this case. In Fig. |B|a 
we have plotted the lowest energy levels (LEL) in each 
sectors Qtot = and Qtot = 7r as a function of W/t in 





mU= 





<^U=- 


-2t 




-4t 




-6t 




-81 




W/t 



FIG. 4. Groundstate energy (E G s) shifted by \N e \U\ 
(a) and real space pairing probability P rs (b) of the 
L = 10 chain at half-filling, as a function of W/t, for 
U/t = 0, —2, —4, —6, —8. The corresponding symbols are indi- 
cated in the figure. The dot-dashed line in (a) represents the 
groundstate energy of the spin-i XY model corresponding to 
the limiting case |W|/t, \U\/t — > oo. 
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FIG. 5. (a) Lowest energy vs. W/t calculated in the 
Qtot = (open squares) and Qtot = 7f (open circles) subspaces 
in the case of a half-filled L — 10 chain at U/t = —4. The 
insets in (a) show the level crossing and singularity in the gen- 
eralized stiffness x(W) = -d 2 E (W)/dW 2 at W c ^ 0.625t. 
(b) Groundstate energy vs. W/t for a half-filled L = 12 chain 
at U/t = —4. The inset in (b) shows the singularity in x(W). 
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FIG. 6. Correlation functions K,(r = L/2) (i= CDW, 
SDW, SC) vs. W/t for the half-filled PKH chain (L = 10) 
at U — —At. The corresponding symbols are indicated in the 
figure. 



the case of a half-filled L — 10 chain and for U/t = —4. 
Indeed, in this case we observe a well defined transition 
from the Q — to the Q = it sector at W c = 0.625i. 
At the same value we observe (see inset in Fig. ||a) a 
singularity in the behaviour of the generalized stiffness 
(GS) 

X (W) - -d 2 E /dW 2 \ u=const 

due to the presence of a kink in the ground state energy 
Eq(W) at the point where the LEL in sectors Qtot — 
and Qtot — n cross each other. When the number of elec- 
trons is N e = 4n the ground state of the system always 
remains in the sector of momentum space with Qtot — 0. 
In Fig. |5|b we have plotted the ground state energy of 
the half-filled PKH model for L = 12. The insets show 
the anomaly in the behaviour of the generalized stiffness 
x(VF) at the same value W c — 0.625i as in the case of 
L = 10 chain. This clearly indicates the irrelevance of 
the finite-size effects already for L = 10 in the consid- 
ered case U/t — —4. We define the critical value of the 
pair-hopping amplitude corresponding to the transition 
into the ry^-paired state at the point where x exhibits a 
singularity. 

To investigate the nature of ordering in the different 
phases we study the behaviour of the correlation func- 
tions. In Fig. H we plotted the absolute values of the 
spin-spin correlator 



•r)>, 



(55) 



the density-density correlator 

Kcdw{t) = j Y,(~V n (Kn)p(n + r)>, (56) 
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where p(n) = p a and the superconducting correlator 



K S c(r) = - ^2(vU n )Vo(n + r)) 



(57) 



at r = L/2 for different values of W/t for the half-filled 
L = 10 PKH chain at U = —4t. We clearly observe the 
existence of two superconducting phases separated by a 
phase with dominating density-density correlations. In 
contrast to generalized stiffness the finite-size effects in 
correlation functions are sufficiently strong for L = 10 to 
use the large-distance properties of the correlation func- 
tions for the precise determination of the W c . However, 
already in the case of the L = 10 chain the essential phys- 
ical details of the phase diagram are comletely transpar- 
ent. 

In order to clarify the nature of the superconducting 
phases corresponding to W > and W < respec- 
tively, we show in Fig. [?] the superconducting correla- 
tor Ksc{f) (a) and the momentum space pairing prob- 
ability P ms (Q) — (AqAq) (b) as a function of W/t, for 
U/t = -4 for the L — 10 chain at half-filling. 

For W < Ksc( r ) smoothly increases with increas- 
ing |W|/t. The distribution P ms (Q) has a peak at 
Q = 0. The weight of this peak continuously increases 
with |W|/t — > oo. The probability to find a 7r-pair is 
almost zero. The inset in Fig. ^a shows Ksc{ r ) at 
W/t = -2.6 (open symbols). At W/t = -2.6 the su- 
perconducting correlations are well established and show 
weak decay with distance. In the opposite case, W > 0, 
Ksc{ r ) and P ms (Q — 0) continuously decrease. The su- 
perconducting correlations are completely suppressed at 
OAt < W < W c ~ 0.75*. The expectation value to find 
a Cooper pair with center-of-mass momentum Q = be- 
comes almost zero at W > W c . However at W > W c 
we observe a drastic change of Ksc(r) and P ms (Q)- The 
superconducting correlations strongly increase and the 
strong peak at Q = ir appears spontaneously in the distri- 
bution P ms (Q) . The probality to find a Cooper pair with 
center-of-mass momentum Q — ir quickly approaches its 
limiting value corresponding to the case W/t — > oo. The 
inset in Fig. @a shows K sc {r) at W/t = 2.6 (filled sym- 
bols). At W/t = 2.6 the superconducting correlations 
are slightly stronger than at W/t = —2.6 and clearly 
show alternating behaviour with periodicity of two lat- 
tice spacing. The data presented in Fig. clearly show 
the tItt ordering at W > W c . 



B. The non- half- filled band case 





W/t 

FIG. 7. The superconducting correlator Ksc{r) (a) and 
the distribution P ms {Q) = (A^Aq) (b) vs. W/t for the 
half-filled PKH chain (L = 10) at U/t = -4. The inset in (a) 
shows Ksc(r) at W/t = —2.6 (open symbols) and W/t = 2.6 
(filled symbols). The dashed and dashed-dotted lines in (b) 
correspond to maximum values approached by the distribu- 
tion P ma (Q) for Q = and Q = ir in the limiting case |W]/t, 
\U\/t ^oo. 



In Fig. g we have plotted P rs as a function of W/t 
for various values of the parameter U/t calculated for the 
L = 12 chain and two particular band fillings v — 1/4 and 
v = 1 /3. As we see the transition into the ^-paired state 
occurs also away from half-filling. Moreover, in contrast 
to the half-filled case, the transition is now very sharp. 
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FIG. 8. Real space pairing probability P rs vs. W/t and 
various values of the parameter U /t for L — 12, N e = 6 (open 
symbols) and N e — 8 (filled symbols). 



To study the phase diagram at v ^ 1/2 we investigate 
the behaviour of the CDW and SC correlations. In Fig. 
||a we have plotted the correlators Kcdw and K$c a t 
r = L/2 vs. W/t calculated for the L = 12 PKH chain in 
the case of v = 1/3 and various values of the parameter 
U/t. As it follows from the behaviour of the correlation 
functions the sector with dominating density-density in- 
stability shrinks to a narrow strip with increasing on-site 
attraction. If for U/t = —4 the density-density corre- 
lations still dominate for W ~ 0.4i, for U/t — —6 only 
superconducting correlations remain. 

In Fig. ^b we show the same correlators calculated 
for a fixed on-site attraction U/t = —4 and two different 
band fillings v\ — 1/4 and i>2 = 1/3. As it is clearly 
seen the decreasing band filling leads to a reduction of 
the range of coupling constants where the system shows 
non-superconducting behaviour. In agreement with the 
strong-coupling expansion analysis of Sect. II we con- 
clude that the system exhibits a transition from the BCS 
type singlet superconducting phase at W < W c into an 
tin paired state at W > W c . 




-3 -2 




w/t 




C. Nature of the transition 

In order to show the finite bandwidth nature of the 
transition into an ^-superconducting state in Fig. [h] 
we have plotted the momentum distribution 

a 

in the ground state of the L — 10 half-filled PKH chain 
as a function of the parameter W/t at U/t = —4. In the 
case of an attractive Hubbard model (W = 0, U = —At) 
the momentum distribution in the groundstate has a 



FIG. 9. Correlation functions KcDw(r = L/2) (filled sym- 
bols) and Ksc(r = L/2) (open symbols) for the L — 12 chain, 
v = 1/4 and U/t = -4, —6, —8 (a) and for the L = 12 chain, 
U/t = —4 and v = 1/4 (squares) and v = 1/3 (circles) (b). 
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FIG. 10. (rife) distribution for different values of the pa- 
rameter W/t, in the case of half-filled L = 10 chain and for 
U/t — —4. The corresponding symbols are indicated in the 
figure. 

shape of the standard Fermi distribution: the states with 
\k\ < kp = tt/2 are almost completely occupied, while 
the states with \k\ > kp = tt/2 are almost empty. In 
the case of W < the momentum distribution remains 
qualitatively unchanged. With increasing \W\/t the dis- 
tribution (nfc) in a monotonic way approaches its limiting 
behaviour (n*,) = 1 at |W|/t — * oo. 

At W > the standard shape in the momentum 
distribution remains for < W < W c \ however, at 
W ~ W c ~ 0.625i we observe an abrupt change in the 
momentum distribution. For W > W c (rife) ~ 1 for all 
k. States with different momenta are almost equally oc- 
cupied by electrons. There is no trace of the Fermi dis- 
tribution. As far as states with momentum k and n — k 
are almost equally occupied after the transition, the con- 
tribution of the one- electron band to the groundstate en- 
ergy is almost completely suppressed. The ground state 
energy of the system becomes linear in W. Therefore we 
conclude that the transition into an %-superconducting 
state corresponds to destruction of the single- electron 
conduction band and the creation of a strongly- correlated 
Tin -pair band. 

As far as the transition into an ry^-paired state is con- 
nected with the formation of a ry^-pair (doublon) band, 
after the transition the characteristic length scale in the 
system £ ~ . Therefore the finite-size effects are very 
weak. In Fig. [ll] we have plotted W c vs. U calculated 
for the half-filled chain with L = 6, 8, 10, 12. As it can be 
seen in this figure the finite size effects have only a very 
weak influence on the critical value. 



FIG. 11. W c vs. U calculated for the half-filled PKH chain 
with L = 6, 8, 10, 12. 

V. SUMMARY 

To summarize, in this paper we have studied the 
groundstate phase diagram of the one-dimensional 
Penson-Kolb-Hubbard model using the continuum-limit 
field theory approach and finite system numerical stud- 
ies. We have shown that the model has a very rich 
phase diagram including, at half-filling (see Fig.||), two 
different superconducting phases and four different in- 
sulating phases corresponding to the singlet supercon- 
ducting phase (U, W < 0) with BCS type pairing at 
|E7|,|W| <§C t, monotonically evolving into the local- 
pair 77 -superconducting phase with U/t and/or W/t ap- 
proaching — oo; the T^-superconducting phase (W > 
W C {U) > 0); the Mott antiferromagnet [U > 2\W\) 
phase, the Peierls dimerized insulator (0 < U < —2W); 
the CDW insulator (W > 0, U < -2W); the unconven- 
tional insulating phase characterized with coexistence of 
the CDW and bond-located staggered magnetization Bd- 
SDW (W > 0,U < 2\W\). 

Away from half-filling (see Fig.[|) the phase diagram 
consists of the singlet-superconducting phase (W < 
min{— U/2, W C (U)}), and the metallic Luttinger-liquid 
(LL) phase {-U/2 < W < W C (U)), and the r? w - 
superconducting phase (W > W C (U) > 0). With increas- 
ing on-site attraction the LL phase shrinks and at a criti- 
cal value U c , determined by the condition U C +2W C (U C ) — 
0, only the critical line W C (U) remains for U < U c . 
In this range of couplings the transition from an 770- 
superconductor to an ri^-superconductor is realized. The 
critical value W C (U) weakly depends on the band filling 
displaying a non- monotonic behaviour: at \U\/t <C 1 a 
slight decrease while at \U\/t^> 1 a slight increase with 
increasing v. 

The obtained phase diagram is in agreement with re- 
sults of previous studies based on the real space renormal- 
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ization group method pq | , the continuum-limit approach 
p6| , the slave-boson mean- field method |3^] as well as for 
the particular case of U = with the results of numer- 
ical studies using exact (Lanczos) diagonalizations |27| l 
and DMRG studies @. 

Using Lanczos diagonalization we have studied in de- 
tail the transition into the ry^-paired state. After the 
transition the on-site pairing probability sharply ap- 
proaches its limiting value corresponding to the case 
W/t — > oo and the ^-superconducting correlations dom- 
inate in the system. The transition corresponds to 
an abrupt change in the groundstate structure. Af- 
ter the transition the single-electron conduction band is 
completely destroyed and the strongly-correlated ry w -pair 
(doublon) band is established. The transition occurs at 
any band filling, the critical value W C (U) weakly depends 
on the band filling but is strongly renormalized by the on- 
site Hubbard interaction: W C (U) is linear in U for \U\/t 
of the order of unity and is inversely proportional to \U\ 
i.e W C (U) ~ -t 2 /\U\ for \U\/t > 1. We believe that this 
phase diagram is generic for the Penson-Kolb-Hubbard 
model and will remain unchanged in higher dimensions. 
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